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RESUMEN 

Se efectuaron tres simulaciones auto-consistentes de A-cuerpos de fusiones 
entre galaxias, con un perfil pronunciado de materia oscura, para estudiar si la 
pendiente interna del perfil de densidad oscura se preserva en los remanentes. En 
estas simulaciones los progenitores incluyen un disco estelar y un cspfn intrmsico 
en el halo a diferencia de estudios similares. La razon de masas de las galaxias 
progenitoras es de alrededor de 1:1, 1:3 and 1:10. Encontramos que el perfil de 
densidad inicial pronunciado de los halos oscuros se preserva en los remanentes 
para los casos considerados aquf. 

ABSTRACT 

We carried out three self-consistent A-body simulations of galaxy mergers, 
with a cuspy dark matter profile, in order to study if the inner dark density slope 
is preserved in the remnants. In these simulations the progenitors include both 
a stellar disk and an intrinsic angular momentum for the halos, unlike previous 
similar studies. The mass-ratios of progenitor galaxies are about 1:1, 1:3 and 1:10. 
We find that the initial cuspy density profile of the dark halos is preserved in the 
remnants for the cases considered here. 

Key Words: GALAXIES: KINEMATICS AND DYNAMICS GALAX- 
IES: DARK MATTER METHODS: AT-BODY SIMULA- 
TIONS. 



1. INTRODUCTION 

High-resolution simulations of hierarchical struc- 
ture formation within the ACDM cosmogony show 
that the inner regions of the halos follow a power- 
law cusp, p oc ?"~ 7 , with 1 < 7 < 1.5 while in their 
outer parts the profiles go as p oc r -3 over a wide 
range of mass-scales (e.g., Navarro, Frenk &; White 
1997, hereafter NFW, Moore et al. 1999, hereafter 
M99, Fukushige, Kawai & Makino 2004, Navarro 
et al. 2004). 

The physical origin of this profile is still unclear 
and it is the subject of an intense investigation. Nev- 
ertheless, there is a certain agreement in the sense 
that such density profile is linked to the accretion 
and merging history of the dark matter substructures 
(e.g., Syer & White 1998, Dekel, Devor & Hetzroni 
2003, Williams, Babul & Dalcanton 2004). 

Several authors have studied the shape and sur- 
vival of the density profile inner slope in collisionless 
mergers (e.g., White 1978, Fulton & Barnes 2001, 
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Boylan-Kolchin & Ma 2004). In the early study of 
White (1978) it was found that mergers of cored 
galaxies were more concentrated than their progeni- 
tors. Fulton & Barnes (2001) found in their A-body 
experiments that steep cusps are preserved during 
mergers of equal-mass progenitors. Boylan-Kolchin 
& Ma (2004, hereafter BKM) concluded that if two 
galaxies merge having each one a core-type profile 
the resulting remnant will be core-type, and if one 
progenitor has a cuspy profile then the remnant will 
also be cuspy. A result in agreement with other sim- 
ulations (Moore et al. 2004). 

These type of studies find justification in that al- 
though in a collisionless simulation the phase-space 
density must be conserved, due to Liouville's theo- 
rem, there is no a priori reason to expect that both 
configuration and momentum distributions are to be 
conserved independently (BKM). 

A common feature of previous related works is 
that they consider gravitational systems with only 
a spherical or dark matter component, and they do 
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not consider the effect, for example, that a disk or 
an intrinsic halo spin might have on the preservation 
of cusps during a merger process. This needs to be 
addressed since there are several studies that indi- 
cate that angular momentum may play a key role at 
the time of formation of the halo as well as on the 
shape of the density profile within the inner regions, 
by preventing that dark matter particles reach the 
inner regions, resulting in a shallower central profile 
(e.g., White & Zaritsky 1992, Hiotelis 2002, Ascasi- 
bar et al. 2004). Furthermore, cosmological simula- 
tions lead to dark halos having an intrinsic angular 
momentum (Barnes & Efstathiou 1987, Lemson & 
Kauffmann 1999). 

In this work we study three numerical A^-body 
simulations of progenitor self-consistent disk galaxies 
with mass-ratios of about 1:1, 1:3 and 1:10, follow- 
ing parabolic encounters. The present work expands 
on earlier similar studies mainly in that the progen- 
itors have a disk component and their halos have an 
intrinsic angular momentum. This would allows us 
to test, using more realistic initial conditions, if the 
cuspy nature of the progenitor galaxies is preserved 
in mergers; at the resolution provided by our A^-body 
simulations. 

The rest of this paper has been organized as fol- 
lows. In §2 we describe the theoretical method used 
to construct the galaxy models, the properties of 
their dark halos, especially the behavior of the in- 
ner logarithmic derivative and the initial conditions 
for the encounters. In §3 we present our results, and 
in §4 we summarize our main results. 

2. MERGER MODEL 

2.1. Disk Galaxies 

Our galaxy models consist of a spherical halo and 
a stellar disk; no bulge component is included. The 
disk is represented by an exponential profile in the 
radial direction and by an isothermal sheet in the 
vertical one (Hernquist 1993): 



Pd{R, z) 



exp(-i?/i? d )sech 2 (z/z rf ) , (1) 



where Rd and z& are the radial and vertical scale- 
lengths of the disk, respectively. 

For the dark halo, a spherical model is adopted by 
assuming a NFW-profile with an exponential cutoff: 



Ph(r) 



exp 



- - 

Woo 



(2) 



where r s and r2oo are the the scale and "virial" radii, 
respectively, c = q~ 1 = r2oa/r s and Mh are the con- 



centration and the mass of the halo. The r2oo is de- 
fined in cosmological simulations as the radius where 
the mean interior density is 200 times the critical 
density. In general, r2oo is not exactly the virial ra- 
dius (§2.4), but is a good approximation to it and 
serves to characterize the halo (NFW, Cole & Lacey 
1996). 2 Here, ah is a normalization constant given 
by 



exp(<7 2 ) 



V^exp(g 2 )Erfc(g) + \ cxp(g 2 )Ei((? 2 ) - 1 



being Erfc(a;) the complimentary error function and 
Ei(x) the exponential integral. 

The radial scale length of the disk is obtained 
accordingly to the Mo, Mao & White (1998) frame- 
work of disk galaxy formation. To determine this 
quantity five parameters are required, namely: the 
circular velocity V c at r2oo, the dimcnsionlcss spin 
parameter A, the concentration c of the dark halo, 
the fraction of disk to halo mass arid, and the frac- 
tion of angular momentum in the disk to that in the 
halo jd- We followed the method described by Shen, 
Mo & Shu (2002) to obtain these parameters. A 
ACDM cosmology is adopted with a mass density of 
fio = 0.3 and with a contribution due to the cosmo- 
logical constant of £!a = 0.7. 

The selected galaxy disks satisfy a stability cri- 
terion e m > 0.9; where e m =V m (GMd/ Rd) -1 ' 2 , being 
V m the maximum rotation velocity (Efsthatiou, Lake 
& Negroponte 1982, Syer, Mao & Mo 1997). Mass- 
ratios of about 1:1, 1:3 and 1:10 were chosen to study 
the properties of the merger remnants. Finally, ve- 
locities of particles were set up through Jeans equa- 
tions following the method prescribed by Hernquist 
(1993). 

In Table ^ we list the parameters defining our 
galaxy models; below the merger simulation label 
(51, 52, and 53) the mass ratio of the secondary to 
the primary galaxy has been indicated. The number 
of particles in the halo A^ and in the disk are 
also listed. 

2.2. Encounter Parameters 

Only parabolic encounters are considered for the 
binary mergers. The relative initial separation be- 
tween both galaxy centers is taken to be 



Rs = 1.25 (r 2 oo,i + ^200,2); 



(3) 



where r2oo,i and r 2 oo ; 2 are, respectively, their r 2 oo 
"virial" radii. 



2 The virial radii of the disk numerical models in §2.4 
turnout to be Ri 20% smaller than the T2oo m Table 1. 
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TABLE 1 
INITIAL GALAXIES 







Halo Properties 






Disc Properties 




MJM V 

St / [J 


M h 


?"200 


A 


3 


c 


A h 


M d 






N d 




[M ] 


[kpc] 










[M ] 


[kpc] 


[kpc] 




51 


4.74 x 10 11 


110.3 


0.112 


0.106 


11.39 


240000 


4.16 x 10 10 


6.6 


1.29 


60000 


0.14 


6.99 x 10 10 


58.3 


0.071 


0.035 


9.94 


165992 


2.88 x 10 9 


1.6 


0.19 


48689 


52 


4.02 x 10 11 


104.3 


0.063 


0.041 


7.63 


57126 


2.11 x 10 10 


2.4 


0.39 


12000 


0.32 


1.25 x 10 12 


152.3 


0.056 


0.068 


3.84 


177571 


7.78 x 10 10 


5.7 


0.47 


44267 


53 


8.11 x 10 10 


61.2 


0.099 


0.142 


10.87 


150000 


7.22 x 10 9 


4.6 


0.73 


30000 


0.98 


8.33 x 10 10 


61.8 


0.122 


0.095 


10.01 


154050 


7.15 x 10 9 


3.9 


0.56 


29686 



Orbital angular momentum is introduced in the 
simulations by randomly choosing a pericenter for 
each encounter, assuming that galaxies are point par- 
ticles in a Keplerian orbit. Pericentric radii in the 
range R p = {5 — 20} kpc were adopted. These val- 
ues for the orbital parameters are consistent with 
those found in cosmological A-body simulations and 
tend to favor mergers (Navarro, Frenk & White 
1995). The relative orientation of the spin vector 
of the galaxies, relative to the orbital plane of the 
encounter, is chosen randomly. 

2.3. Numerical Tools 

These A-body simulations were performed us- 
ing a parallel version of GADGET, a tree base 
code with individual timesteps (Springel, Yoshida 
& White 2001). The runs were made in a Pen- 
tium based cluster of 32 processors (Velazquez & 
Aguilar 2002). Softening parameters £d = 35 pc 
and eh = 350 pc for disk and halo particles, respec- 
tively, were used. Since GADGET uses a spline 
kernel for the softening, the gravitational interaction 
between two particles is fully Newtonian for sepa- 
rations larger than twice the softening parameters 
l|Power et al.ll2003|) . 

The typical time of arrival of the galaxies to the 
first passage through pericenter is about 1 Gyr. We 
follow each encounter for a total time of about 8 Gyr. 
At this time the remnants had reached a stable value 
close to virial equilibrium. Energy conservation was 
better than 0.25% in all our simulations. 

2.4. Resolution Criteria and Stability of Initial 
Dark Halo Profiles 

The two-body relaxation time-scale imposes an 
inner 'convergence' radius r c over which the stellar 
system can be adequately described by a collision- 
less distribution function (e.g., Power et al. 2003; 
Hayashi et al. 2004; Diemand et al. 2004; Bin- 
ney 2004). 



Here, we take r c as the radius where its local two- 
body relaxation time-scale t T is equal to the period 
T v of a circular orbit at the virial radius; i.e., r c 
satisfies 



A 



*r(r) 

T v 8\nNV(r)r v 



= 1 



(4) 



where V v the circular velocity at the virial radius, r v , 
and N — N(r) the total number of particles inside a 
radius r. The virial radius is r v = GM t /\W\, being 
M t the total mass of all bound particles and W its 
total potential energy. 

We computed spherically averaged density pro- 
files, p(r), for each initial galaxy in a logarithmi- 
cally spaced grid within r c and r v ; in general, we 
find r c /r v rj 0.01. To check that the dark den- 
sity profile of our initial galaxy models does not 
change, we evolved in isolation each galaxy for about 
8 Gyr, that corresponds to ss 25 dynamical times, 
= -\/37r/16G/5h; where ph is the mean density 



id 



at the half-mass radius. At this time the virial ratio 
fluctuates within < 1% from the ideal value of unity, 
and no significant change in the density profiles was 
observed; see Figure 1. 

Results in the literature (e.g., Kazantzidis, Ma- 
gorrian & Moore 2004, Springel et al. 2005) indi- 
cate that obtaining the velocities of particles us- 
ing the Hernquist (1993) procedure, that assumes a 
Maxwellian local velocity distribution, is inadequate. 
In these works it is found that models evolved in iso- 
lation tend to relax rapidly to a shallower inner den- 
sity slope than the initial condition. However, our 
models constructed using Hernquist 's method do not 
show a shallowing of the inner slope when evolved in 
isolation (see Figure 1). 

We fitted NFW density profiles, 



P*(r) 



Ps 



(r/f s )(l + r/f s y 



(5) 
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Fig. 1. Density profiles of the dark matter component 
at different times of a galaxy model evolved in isolation. 
Different lines indicate the time spent in isolation. The 
dynamical time is idyn ~ 0.3 Gyr, so the evolution is for 
~ 25tdyn- The inner profile, close to r c , remains stable 
showing no shallowing of the density slope. The inset 
indicates the evolution of the virial ratio. 

to the p(r) of the initial numerical models; where 
f s is the scale radius. The fitting was done 
by Y 2 -mi nimization using t he Levenberg-Marquardt 
method l|Press et alJ ^92) . From the fitted pro- 
file we computed the logarithmic derivative (3{r) = 
dlog p/dlogr. 

To verify that our binning method does not affect 
the fitting of the density profiles, we computed p(r) 
in bins of equal number of particles, from 100 to 
1500, and the parameters of the fitted profile did 
not changed significantly (< 1%) in comparison to 
the logarithmically spaced grid. We consider this 
provides a good degree of confidence of the results 
presented for our fits. 

The mean RMS of the deviations, S = (p — p*) / p, 
for the initial profiles in the complete fitting interval 
is 0.120, while for the inner region up to a radius 
of 10% the virial radius a mean, RMSio, a value of 
0.078 is obtained. The latter value is an indication 
that the iV-body realization follows closely a NFW 
profile at the inner region. 

We estimated the logarithmic derivative at r c , 
/3 C , for each individual galaxy, and a mean value of 
— 1.21 is found. Individual values of (i c and RMS, 
both total and inside 0.10i? v , are listed in Tabled 
We note that the values of (3 C differ from each other 



TABLE 2 
INITIAL GALAXY SLOPES 
Merger TV-body Model 



M s /M p 


RMS 


RMSio 


-ft 


51 


0.098 


0.046 


1.25 




0.115 


0.068 


1.22 


52 


0.128 


0.105 


1.28 




0.178 


0.120 


1.06 


53 


0.102 


0.075 


1.12 




0.100 


0.053 


1.32 



since r c for each progenitor is different. 

3. MERGER DENSITY PROFILES 

In Figure El we show the time evolution of the 
dark matter density profile for merger 52, from 
6 Gyr to ~ 8 Gyr. As observed, the inner density 
profile does not change significantly from sa 6 Gyr 
to the end of the simulation. The evolution of the 
virial ratio for the whole time of the simulation is 
shown in the inset. The dynamical time of the rem- 
nant is tdyn ~ 0.4 Gyr, so for ss 5t dyn the profile re- 
mains unchanged and the virial ratio has settled to 
the equilibrium value around 1; with some numerical 
fluctuations less than 1%. A similar behavior in the 
profile and virial ratio was observed for the other two 
simulations considered. The analysis of the resulting 
profile of the remnants were done at t 8 Gyr. 

A NFW profile was fitted to the resulting merger 
remnants. The fitting is done in their corresponding 
interval (r c ,r v ). In Table 01 we present the values 
for the total RMS deviation of the fit (column 2), 
the RMS inside 0.1r v (column 3), the scale radius 
f s in kpc (column 4), the concentration parameter 
c = r v /f s (column 5), and the logarithmic derivative 
at r c obtained from the fitting (column 6). 

In Figure|31(fop) we show the fractional residuals 
of the merger fittings to a NFW profile and (bottom) 
the numerical logarithmic derivatives of the data 
that correspond to the fitted profile JSJ. Though the 
numerical logarithmic derivatives are rather noisy, 
especially at the boundaries of the fitting region, 
they follow in average the (3(r) obtained from the 
fits. The mean RMS of the density fits for the entire 
interval is 0.093, while for the region from r c up to 
0.10r v is about 0.066. On other hand, the average 
value of p c obtained for the remnants is —1.24. 

The mean RMSio m the fits of the remnants is 
about the same as that found in the fit for the pro- 
genitors to the NFW profile. The values of /3 C of the 
remnants (Table 3) do not necessarily need to coin- 
cide with, say, the average value of their progenitors 
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Fig. 2. Late time evolution of the density profile of the 
remnant of simulation 52, from « 6 Gyr to ~ 8 Gyr. The 
dynamical time of the remnant is td yn ~ 0.4 Gyr. The 
inset shows the virial ratio during the total time of the 
simulation. 

TABLE 3 
REMNANTS' NFW PROFILE. 



Merger 




Profile Parameters 






RMS 


RMSio 


r s 


c 




51 


0.074 


0.040 


5.69 


13.17 


1.24 


52 


0.124 


0.078 


24.40 


6.39 


1.12 


53 


0.082 


0.081 


4.42 


14.96 


1.37 



(Table 2) since r c is different in both cases. We con- 
sider the RMS io of these fits to be a better indicator 
of the preservation of the initial cuspy profile at the 
inner region. 

There has been some discussion (e.g., Navarro et 
al. 2004) if a profile of the form proposed by M99, 



p(r) 



Pm 



(r/r M ) 1 - 5 [l + (r/r M ) 



1.51 



(6) 



provides a better fit to the dark halos. We tested 
also if © provided a better fit to the density pro- 
files of the remnants. In Tabled we list the values 
obtained for similar quantities as shown in Table |31 
but with c = r v /rM- In Figure 0] we show the frac- 
tional residuals of the fits and /3(r) using ©. 

We find a mean f3 c = —1.53 but with the av- 
erage RMSio to be 0.137. A value that is w 100% 
larger than that obtained with the NFW profile in- 
side 0.1i? v - We take this result as an indication that 



0.4 
0.2 

o 

-0.2 
-0.4 

-1 



oa. -2 



-3 



Mergers 





0.01 



0.1 
r/r v 



Fig. 3. (Top) Fractional residuals, <5nfw = (p— /5nfw)/p, 
of the mergers' density profiles with respect a NFW pro- 
file. (Bottom) Numerical logarithmic derivatives (dotted 
lines) and the ones obtained from the fitted profile (con- 
tinuous lines). 

TABLE 4 
REMNANTS' M99 PROFILE. 



Merger 




Profile Parameters 






RMS 


RMSio 




-/3c 


51 


0.084 


0.093 


10.56 7.09 


1.53 


52 


0.196 


0.219 


47.03 3.32 


1.51 


53 


0.077 


0.099 


7.98 8.29 


1.56 



the M99 profile does not provides a good fit to our 
merger remnants. 

4. SUMMARY 

We have built up disk galaxies using the cosmo- 
logically motivated model of MMW to study the cus- 
pyness properties of merger remnants. The progeni- 
tor galaxies have spin in their halos and mass-ratios 
of about 1:10, 1:3, and 1:1. This work was restricted 
to study binary mergers following parabolic encoun- 
ters. 

Our results indicate that angular momentum, 
both intrinsic and orbital, do not change the behav- 
ior of the density profile inside the inner regions of a 
merger remnant. 

However, taking into account that few simula- 
tions were done, we can not completely rule out the 
effect of the initial angular momentum in the sys- 
tem on the inner density profile of the remnants. A 
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Fig. 4. (Top) Fractional residuals, 8m99 = (p — PM99)/p, 
of the mergers' density profiles with respect a M99 pro- 
file. (Bottom) numerical logarithmic derivatives (dotted 
fines) and the ones obtained from the assumed profile 
(continuous lines). 

more wider set of simulations is required to check the 
result found here. 

Acknowledgments 

We thank CONACyT, Mexico, for financial sup- 
port through Project 37506-E. An anonymous ref- 
eree is thanked for providing important comments 
that help to improve the presentation of this work. 

REFERENCES 

Ascasibar Y., Yepes G., Gottlober, Miiller V., 2004, MN- 

RAS, 352, 1109 
Barnes J., Efstathiou G., 1987, ApJ, 319, 575 
Binney J., 2004, MNRAS, 350, 939 
Boylan-Kolchin M., Ma C, 2004, MNRAS, 349, 1117 

(BKM) 

Cole S., Lacey C, 1996, MNRAS, 281, 716 
Dekel A., Devor J., Hetzroni G., 2003, 341, 326 
Diemand J., Moore B., Stadel J., Kazantzidis S., 2004, 

MNRAS, 348, 977 
Efstathiou G., Lake G., Negroponte J., 1982, MNRAS, 

199, 1069 

Fukushige T., Kawai A., Makino J., 2004, ApJ, 606, 625 
Fulton E., Barnes J. E., 2001, Ap&SS, 276, 851 
Hayashi E., Navarro J. F., Power C, Jenkins A., Frenk 

C. S., White S. D. M., Springel V., Stadel J., Quinn 

T. R., 2004, MNRAS, 355, 794 
Hernquist L., 1993, ApJS, 86, 389 
Hiotelis N., 2002, A&A, 382, 84 

Kazantzidis S., Magorrian J., Moore B., 2004, ApJ, 601, 
37 



Lemson G., Kauffmann G., 1999, MNRAS, 302, 111 
Mo H. J., Mao S., White S. D. M., 1998, MNRAS, 295, 
319 (MMW) 

Moore B., Quinn T., Governato F., Stadel J., Lake G., 
1999, MNRAS, 310, 1147 (M99) 

Moore B., Kazantzidis S., Diemand J., Stadel J., 2004, 
MNRAS, 354, 522 

Navarro J. F., Frenk C. S., White S. D. M., 1995, MN- 
RAS, 275, 56 

Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 
490, 493 (NFW) 

Navarro J. F., Hayashi E., Power C, Jenkins A. R., Frenk 
C. S., White S. D. M., Springel V., Stadel J., Quinn 
T. R., 2004, MNRAS, 349, 1051 

Press W. H., Teukolsy S. A., Vetterling W. T., Flannery 
B. P., 1992, Numerical Recipes: The Art of Scientific 
Computing. Cambridge Univ. Press, New York 

Power C, Navarro J. F., Jenkins A., Frenk C. S., White 
S. D. M., Springel V., Stadel J., Quinn T., 2003, MN- 
RAS, 338, 14 

Shen S., Mo H. J., Shu C, 2002, MNRAS, 331, 251 

Springel V., Yoshida N., White S. D. M., 2001, New As- 
tronomy, 6, 79 

Springel V., Di Matteo, T., Hernquist L., 2005, MNRAS, 
361, 776 

Syer D., White S. D. M., 1998, MNRAS, 293, 337 
Syer D., Mao S., Mo H. J., 1997, astro-ph/9711160 
Velazquez H., Aguilar L. 2003, RMxAA, 39, 197 
White S. D. M., 1978, MNRAS, 184, 185 
White S. D. M., Zaritsky D., 1992, ApJ, 394, 1 
Williams L. L. R., Babul A., Dalcanton J. J., 2004, ApJ, 
604, 18 



G 



ACEVES & VELAZQUEZ 



Hector Aceves & Hector Velazquez: Institute de Astronorma, UNAM Campus Ensenada, Apdo Postal 877, 
Enscnada 22800, Baja California, Mexico (aceves@astrosen.unam.mx, hmv@astrosen.unam.mx) 



